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ABSTRACT 

This paper brings together a number of algorithms for the estimation of the envelope and 
instantaneous frequency functions suitable for machinery diagnostics and structural health 
monitoring. This includes the classical approach based on Rice's frequency, the most 
commonly used method based on the Hilbert transform and recent developments related to 
the Wigner-Ville distribution and wavelet analysis. The study includes two application 
examples from fault detection in gearboxes and identification systems. The paper shows that 
the envelope and instantaneous frequency functions are an Important link between classical 
Fourier approach based methods and time-variant analysis. 
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1. INTRODUCTION 

The viewpoint of this paper is that there exist three major areas within machinery 
diagnostics and structural health monitoring where the instantenous characteristics, namely 
the envelope and instantaneous frequency, are used. These are: damage detection in rotating 
machinery, identification of nonlinear systems and deterministic approach to stationarity. 

Machines and structures in operation generate forces and motions that produce 
different types of vibration. Any fault in machine operation or structural damage can be 
considered as an additional excitation which results in a change of vibration response. Often 
the overall process of vibration response exhibit modulations, intermodulations or nonlinear 
behaviour due to operational nature or fault. Identification of modulation sources can provide 
valuable information of many faults in rotating machinery, particularly in gearboxes and ball- 
bearings. Modulations in gearboxes are mainly caused by the tooth meshing process, due to 
varying stiffness, or by the eccentricity and local tooth faults 1 . This is exhibited by modulation 
sidebands in vibration spectra, 

Intermodulations results from sum and difference frequencies of low frequency 
harmonics of the shaft speed, harmonics of the tooth meshing frequency and modulation 
sidebands. Modulations in ball-bearings are mainly caused by locaJ surface defects which 
produce impulses often leading to resonances. For advanced defects, simple parameters 
such as kurtosis or crest factor 2,3 are sufficient to detect damage. Often modulation and 
intermodulations which result from damage can be identified by means of ZOOM-spectrum 
and cepstrum 4,5 . For complex machines the signal demodulation procedure, or in other words, 
the analysis of the envelope and instantaneous frequency functions, is very effective to 
identify local tooth damage 6 " 10 and local defects in ball-bearings 1 1 " 14,9 " ig . The amplitude 
demodulation procedure is often called envelope detection or, when used in ball-bearings, the 
high-frequency resonance technique. 

The second important area where the envelope and instantaneous frequency 
functions are often used is identification of nonlinear systems. Since nonlinearities may result 
from damage in a system, the distinction between linear and nonlinear behaviour is an 
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important problem in damage identification. Many different procedures for the identification of 
nonlinear systems from vibration test data have been developed in structural dynamics. It is 
well known that many types of nonlinearities cause a varying nature for the restoring forces 
and natural frequencies of the system. This varying nature of vibration can be studied using 
the envelope and instantaneous frequency functions which lead to nonlinear back-bone and 
damping curves. 



The overall vibration produced by machinery and structures is either stationary or 
nonstationary. The analysis of the envelope and instantaneous frequency functions often 
forms the deterministic test for stationarity in vibration analysis 15 . In this approach, the 
Q process is said to be stationary if it consists of components such that their instantaneous 

n 

■ssr 

ry behaviour, as described by the envelope and instantaneous frequency, does not depend on 



time. Conversely, the process is non-stationary if its instantaneous characteristics depend on 



'y time. The instantaneous characteristics are used in machinery diagnostics to detect local 

a faults. This fault detection method is in fact a test of stationarity and any departure of 

fy instantaneous characteristics from time independent functions is considered as a symptom of 

yf a fault in the system. 

LJ, 

o 

Recent developments in signal processing show that the envelope and instantaneous 
frequency functions are an important link between classical Fourier approach based methods 
and time-variant analysis, as shown in 15 . There exist a number of procedures, developed 
over the last twenty years, which can be used to estimate the envelope and instantaneous 
frequency functions. It appears that the Hilbert transform based procedure 16 is the most 
common approach to the problem. Recent developments include the Wigner-Ville 
distribution 17 and the wavelet transform 18 . However, there exist a number of different 
algorithms which can be used for single- and multi-component signals (frequencies). 

The aim of this paper is to bring together the algorithms of signal demodulation. It is 
hoped that this form of presentation will help to implement the most suitable approach for a 
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given vibration problem. However, the paper does not intend to survey various application 
examples. 



The structure of the paper is as follows. Section 2 shows the complexity of 
modulation and intermodulation processes which can be found in rotating machinery. The 
algorithms of the amplitude and frequency/phase demodulation are given in Sections 3 and 4, 
respectively. This is followed by Section 5 with examples related to fault detection in 
gearboxes and identification of nonlinear systems. Finally, the paper is concluded in Section 

Q 

O 2. MODULATION AND INTERMODULATION PROCESSES IN MACHINE DIAGNOSTICS 

«1 In machinery, modulation is the variation in the value of some parameter which 

characterizes a periodic osciliatjon. This process is an effect of either design features of 

fjj machines (e.g. crankshafts) or more often a result of a fault. Thus the identification of 

q modulation sources can provide valuable information about possible faults in machinery. 

~: Amplitude modulation processes, in general, are caused by impacts and impulses in the 

structure, generated by different phenomena, for example, cracks or missing teeth in gears, 
point defects in ball-bearings, blade cracks in fans or any other local faults. The impulses are 
repeated periodically for each revolution of the wheel, shaft, etc. The frequency modulation 
processes are an effect of parametrical excitation in the structure (e.g. fluctuation in tooth 
loading in gearboxes, varying stiffness of teeth in gears, varying cross-section of shafts, etc.). 
More details can be found in 1,19 " 25 . In what follows some spectral properties of modulated 

processes will be relevant to machinery vibration as summed. 

> 

The vibration signal generated by a machine can be analytically described as a complex wave 

xif)=jj\ + a r (r )V (<)] + „( t ) 

0) 
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where a r (t) is the slow-varying amplitude modulation process, $ r (t) is the fast varying 
phase modulation process, f or are the carrier frequencies (6> or =2x f of -const) and n{t) 
is the noise. Amplitude and phase modulation processes can be represented by, 



(2) 



O) 



where M ir are the coefficients of amplitude modulation intensity and, /}& are phase 
deviations. Substituting (2) and (3) into Eq. (1), it follows that 



/=1 



i jtocr*+ Z /Ar cos ^ r 

■ *=1 +«(*) 



(4) 



Eq. (2) can be replaced by, 



(5) 



The well known Fourier-series expansion, 



(6) 

where J n {fijcr) is the Bessel function of the first kind, together with Eqs. (4) and (5) yields, 
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) in Eq. (6) is the Besset function of the first kind. Eq. (7) can be written as, 
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(8) 



or finally. 



/Hv +00 , \ 



(9) 



Analysing Eq. (9) one can notice that in the frequency domain the vibration signal 
x(t) is represented by R families of spectral components. Each family yields five types of 
side-frequency terms: 



A i e^ or 1 - carrier frequency 
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A 2 e Mor±&ir)t _ sid6 f requenc ies due to amplitude modulating waves 
A e A<»or^ a kr)t _ sjde frequencies due to phase modulating wave 



j(eo or ineo\ r ±ne)2 r ± .... ±na}^ )t _ beat side f re q Ue ncies due to phase modulation 



. a j{t» <)r ±Cl }r ±Ci.2r±-- ±^ir±»^\r^n(O2 r ±...±n0ia.)t . beat side frequencies due to 
amplitude and phase modulation, 

where ,4, are the Bessel function components according to Eq. (9). Fig. 1 shows an example 
of the amplitude and phase modulated signal, which includes two families of spectral 
components, and its spectrum. 

In the case of the pure amplitude modulation process (av =0)the total energy of the 
signal will increase. The energy of the carrier remains unchanged and the additional energy 
which comes from the modulation wave will be represented in sidebands. This process yields 
only one pair of sidebands: sum and difference components. In comparison the total energy 
of the phase (frequency) modulated wave remains constant. To achieve this the energy of the 
carrier will decrease and the difference will be distributed to its infinite number of sidebands. 
In addition phase (frequency) modulation can produce beat components, which are not 
characteristic for the amplitude modulation. In practice, for example, the carrier frequency in 
the gearbox deals with the meshing frequency. Thus one should remember that greater 
advancement of a fault does not need to cause the increase of the level of the meshing 
frequency. This means that the level of the meshing frequency might not be a good indication 
of the advancement of a fault. Additionally, if one considers the phase angle between the 

amplitude and phase modulation waves i.e. a^s^M^. cos(fl,> r+f2,> 0 ). thus the 

i=l 

asymmetric distribution of sidebands will be observed. These brief considerations show that 
the spectrum of the vibration signal generat d by a complex machine might be very 
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complicated. Fig. 1 shows an example of amplitude and phase modulated signal and its 
spectrum. 

3. AMPLITUDE DEMODULATION ALGORITHMS 

Amplitude demodulation is the process of extracting slow variation of the amplitude 
from the modulated signal. This operation falls into two broad categories of synchronous 
detection and envelope detection 26 . Both categories are summarised below. Additionally, the 
algorithm involving Discrete Fourier Analysis (DFA) 27 is also described. 

3.1. Synchronous Detection 

Synchronous detection uses the operation of frequency conversion. The concept is 
outlined in Fig. 2. Here the oscillator is assumed to be exactly synchronised with the carrier 
frequency. 

Assume an amplitude modulated signal x(i)=a(t)cosa Q t where co 0 is the carrier 
frequency. Multiplying both sides by 2zosG> 0 i one obtains 

2x(r)cQ$a 0 t=a(t)+a(t)co$2& 0 1 

(10) 

The first term on the right side is the one we want, since it is an amplitude modulation 
signal. The second term can be removed by low-pass filtering. This procedure requires the 
carrier frequency determination in order to translate the signal spectrum to a new carrier 
frequency equal to zero. Nevertheless some degree of asynchronism must be expected, 
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since the carrier frequency cannot always be exactly determined (for example in the case of 
varying rotational speed of the gearbox). This means that instead of multiplying by 2cosa> 0 t 
one may use the value 200$^ +a>' 0 )t which results in, 

2 x(t)cos(a 0 +<D' 0 )t - a(/)cos<y; f+a(r)cos(2® £ , +a>' 0 )r 

(11) 

Here the second term can still be removed but the harmonic relationships in the 
modulation signal a{i) are destroyed. Therefore instead of a(t) one receives a(t)cos<»' 0 t , 
which in the frequency domain, according to the modulation theorem 28 , gives us the shifted 

n spectrum j. without any possibility of correction. The second difficulty with this 



N procedure is related to the digital filtering operation which reduces the efficiency of the 

algorithm. 



■ An alternative implementation of this method is presented in section 4.2, where the 

fU frequency modulation signal <j>{f) is assumed not to be equal to zero, which allows one to 

O compute either a(t)or 4>{t) 

H 

3.2. Envelope Detection 

There exist a number of envelope definitions. All of them come from random vibration 
or communication systems theory. Physically, for the narrow-banded process x(t). the 
envelope is a smooth curve joining the peaks of x(t). A few envelope definitions and their 
statistical properties have been described by Langley 29 . In what follows a few possible 
implementations of envelope detection methods are briefly discussed. 



3.2.1. Local maxima distribution 
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The envelope as a local maxima distribution was defined by Crandall as being a 
smooth gradial curve joining the peaks and the expected time spent by the envelope between 
the level a and a+da is just the number of those peaks, whose magnitude lies between a 
and a+da multiplied by the expected period for cycles of amplitude a. In practice the 
method leads to running averages of signal local maxima. If the modulated process is written 
as x(t)=a(t)co$[a 0 t+$(t)], the running averages can be carried out as 32 

. e(t)t*E r {Mac\x(t]}=j f Max{a(t)QOs[a 0 1 + ^{f%dr 

(12) 

where E T is a symbol of averaging. If the phase process Is constant (&(t)- const) , for the 
optimal averaging time z» — . Eq. (12) can be replaced by 

•CO-i' 7*(0<*««(<) 

T w 
T 

(13) 

This definition implies the bandwidth restriction. The efficiency of the method depends 
on the choice of the optimal averaging time. It also increases if the process is more narrow- 
banded. The local maxima finding procedure can be done by the approximate construction of 
tangents structure. Unfortunately the envelope is then given as a series of finite number 
samples, which makes exact tangent points calculations impossible. Additionally, the 

33 

integration procedure requires a few samples in a period to reach a good accuracy 
3.2.2. Energy based distribution 

The energy based distribution envelope proposed by Crandall 30 " 31 is related to 
stochastic processes, if the stochastic equation of motion is assumed to be of the form" 0 
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x+/M(*,x)+G(*)=.F(0 

(14) 

where x is a displacement response process. A is assumed to be an odd function with 
respect to x. ft is a small parameter. F{t) is the white noise excitation, the envelope 
function a{t) of the random process x(t) is given implicitly by 30 



1^2 

(15) 



where Jc(<)^ • The right-hand side of the above equation is the sum of the kinetic 

o 

and potential energy per unit mass. In the linear case one can set fim2&> e , F{x,x)=x , 
G{x)=wl x and Eq. (15) can be replaced by 



a(0= 



x 2 (/) + 



'm 2 



c 



(16) 

where co c is some constant frequency. Langle/ 9 has suggested to use <o c as the mean zero 
crossing frequency given by <o c - ^ . where m n is the n^th spectral moment of the single 
sided spectrum of x(t) defined as 

00 

0 

(17) 

This definition of the envelope does not imply any bandwidth restriction. It is easy to 
notice that such an envelope follows the peaks of x(r). since the envelope and the process 
coincide at a peak, where x{t)=0 . The implementation of the method can be realized 
following Eq. (16). The effect of the algorithm is shown in Fig. 3. This requires signal 

24 
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differentiation, which is not a simple task 33 The differentiation can be avoided by measuring 
acceleration or velocity where possible and integrating to obtain velocity or displacement The 

33 

accuracy still requires sampling at over a few times the highest frequency of interest . 

Langley 29 has shown that the mean rate at which the envelope crosses a given level 
with positive slope depends on the 4-th spectral moment m 4 . He indicated that some random 
processes, such as for instance the response of a linear system to white noise, have a 
theoretically infinite value of m 4 . This predicts that the envelope crossing rate is infinite. In 
practice signals generated by machines have a finite value of m 4 . if not, it is still possible to 
filter the response spectrum at a frequency which yields a finite value of m A . 

3.2.3. Rice's envelope 

The so-called classical definition of the envelope has been given by Rice 34 . Writing 
the modulated signal in the form 

n 

(18) 

and selecting a frequency f m called the "midband frequency" Eq. (18) can be replaced by 
n 

(19) 
where, 

l C =Z C « cos K flJ « -/™) f + *n (01 
« 

(20) 
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The expression 



aif) = V^C + /| 

(21) 

is termed by Rice the envelope of x(t) referred to frequency /„ . It has been shown by 

Dugundji 35 that this envelope is completely independent of the midband frequency fm- The 
explicit calculation of Eq. (21) is rather impossible. Thus the Rice's envelope is presented only 
as a theoretical case. 

3.2.4. Hilbert transform approach 

The modulated signal x{t) can be replaced on the basis of the analytic signal 

*H<tT 

(22) 

where x H (i) is the Hilbert transform of x(t) being 2 



-28 



(23) 

Thus the envelope suggested by Dugundji 35 is possible 



Of 
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(24) 

It has been shown 2935 that this envelope and Rice's envelope are equivalent and the 
proof of equivalence does not depend upon x{t) being Gaussian. It is easy to notice that for 
x(t) being harmonic (*(/)= A cos * 0 t) the energy based envelope (16) and the envelope 
defined by (24) are also equivalent However, the central carrier frequency of Dugundji's 
envelope is independent of the carrier frequency. 



nj 

M 
si 

RJ 



Langley 29 has assessed when the envelope given by Dugundji (24) will follow the 
peaks of x(t) . The mean and variance of x H {t) when x(t) has a specific value is given by 29 



(25) 



var 



= m 0 q 



(26) 



where q is a parameter which measure the extend to which x(t) is narrow-banded 



2 i m \ 
q = 1 — 



m 0 m2 



(27) 



Thus a{t) given by (24) will follow the peaks of x(t) (*«o) if q is small {x H =0) 
which requires x(t) to be the narrow banded process. It has been shown that the mean rate 
at which the envelope (24) crosses a given level with positive slope is always finite, which 
does not take place in the energy based envelope. 
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The signal processing implementation of the envelope proposed by Dugundji can 
realised according to Eq. (24). Th Hilbert transform definition given by Eq. (23) can 
written in the convolution form 28 

H[x(t)]=x H (t) = ±x(t)* 1 - 

(28) 



Taking into account the signum function, 



sgn*=H 



[-1 /<0 
1+1 *>0 



(29) 



w 
H 

Mi 



and the Fourier transform 28 , 



F[x(t)hx(fh \x{f)e^dt 



-CO 



(30) 



the following is obtained 



(31) 



{__L_1 



(32) 



Thus Eqs. (28), (31) and (32) yield 

F[x H (01= X T (/)=*(/)(-/ sgn /) 

(33) 



28 
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where X T is the signal X{f) with shifted phase by | for the negative frequency 
components and for the positive frequency components. It means, that the Hilbert 
transform x H (t) of the vibration signal *(/) oan be easily obtained by calculating the 
complex spectrum of the signa \(x{f)), shifting the phase (x T (f)) according to Eq. (33) 
and employing the inverse Fast Fourier Transform (FFT). The computation procedure can be 
realized also in another version. Going back to Eq. (22) and using the Fourier transform (30) 
one can find 



X(f)=X(f)+jX T (fhX(f)+s^fX(f)= 



o /<o 
x(f) /=o 

2X{f) f>0 



(34) 



flj This equation represents the spectrum of the analytic signal and can be very easy 

U computed using the spectrum of the signal x(t) . It allows the computation of the envelope 

jy function directly from Eq. (22) as a modulus of the analytic signal. 

3 The method presented in this section seems to be very effective and is often given as 

an example of an effective algorithm of the amplitude demodulation 37 - 21 ' 38 9 . It has been used 

37 21 38 9 32 

in a number of diagnostics applications ' 



3.2.5. Discrete Fourier analysis (DFA) algorithm 



This algorithm has been formulated and tested by Hsueh and Bielawa 27 . The 
assumption is that the amplitude and phase are very slowly relative to the carrier frequency, 
which means that signals are narrow-banded. If they are additionally periodic, they can be 
represented by a time-variable Fourier series. Thus a typical amplitude modulated signal can 
be written as follows 27 



"70 
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x(t)=x 0 x mc (*)cos ca 0 1 + x m (Osiu 6) 0 1 + higher harmonic terms + noise 

(35) 



ru 
ru 



where a> 0 is the carrier frequency, x 0 » are the t' r "e- variable Fourier coefficients of 
*(f), representing the envelope functions. According to the discrete time-variable version of 
the Fourier analysis. These harmonic coefficients are given by 27 



1 '* 

^(0=— Jx(r)cos<o 0 /rfr 



*s 

111 (36) 



where T c is the carrier period and t k is the time instant. Thus in order to detect the 
amplitudes of the envelope parts of the input signals, it is necessary at first to remove, by 
means of analogue and/or digital initial filtering, the higher harmonic terms and high frequency 
noise from Eq. (36). Then the DFA analysis according to Eq. (37) can be performed. It is easy 
to notice that this algorithm is similar to the synchronous detection method described in 
section 3.1. Instead of low-pass filtering after frequency conversion, the DFA algorithm has 
been used. The implementation of this method 27 has shown its very good accuracy and timing 
characteristics. However, a large number of samples in one carrier period is also required 
because of a numerical integration. This method also needs the carrier frequency 
determination. Thus some degree of the carrier frequency shift must be expected and the 
analysis similar to that one presented in section 3.1 can be performed. 



4. FREQUENCY DEMODULATION ALGORITHMS 



•3A 
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The frequency demodulation procedure is connected with the computation of an 
instantaneous frequency /(/) of the signal x(i). According to Eq. (1) this operation can be 
written as, 

(37) 

where 6{t) is the total instantaneous phase of the signal and /o is the carrier frequency. 
Since this operation requires signal differentiation, in many cases frequency demodulation is 
H replaced by phase the demodulation procedure which gives 

N (38) 

■3 

n \ where <a 0 s=2*r/ 0 , Eq. (38) shows that phase demodulation requires removing of the 



iU 



constant rotational term 6>§t related to the carrier frequency. Thus the required earner 

frequency can be determined. Both approaches are similar from the diagnostic point of view 
since they give similar information about the fault. In contrast to amplitude demodulation, 
frequency (phase) demodulation is a nonlinear operation because of the nonlinear 
relationship between the signal x(t) and its instantaneous frequency (phase). Thus the 
procedure is much more difficult. 



There exist many different algorithms of frequency and phase demodulation. Almost 
all of them are included in one of the following categories 28 : (a) FM-to-AM conversion, (b) 
phase-shift discrimination, (c) zero-crossing detection and (d) frequency feedbackA signal 
processing implementation of these algorithms are presented below. 

4.1. FM-to-AM Conversion 
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Any device whose output equals the time derivative of the input produces frequency 
modulation-to-amplttude modulation conversion 28 . Consider the frequency modulated signal 

*(/)==4cos[a> 0 f+^(/)] 

(39) 



Then 



=2 J(f 0 + i-^ccs[^ + 4±180 0 )] 



(40) 

Thus taking into account Eq. (38) one can obtain 

(41) 

where Env\x(t)\ is the envelope function of the differentiated signal *(*)- Input si 9 na[ ' s 
assumed to have some constant amplitude A . In practice an amplitude limiter is necessary at 
the input to remove any variations. Such a limiter is not easy to implement The envelope 
function can be calculated by means of any algorithm presented in section 3. Additionally 
signal differentiation is required. This can be avoided, by analogy to section 3.2.2, by 
measuring acceleration and integration to obtain velocity. This method has been applied in 
machinery diagnostics 10 ' 32 . 



4.2 Synchronous Detection 

Synchronous detection used in section 3.1 to compute amplitude modulation signal 
can be, by analogy, used also in the case of phase modulation. Considering 
*(r)=y4cos[^o' + 0(O] and multiplying both sides by 2cosco 0 t one obtains 
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2 *(/)cos a> 0 r a A cos ^(*)+ ^4 cos[ 2*» 0 f + ^(*) ] 

(42) 

The second term can be, by analogy to section 3.1, removed by low-pass filtering. 
Thus the instantaneous phase process is easy to extract The method used simultaneously 
for amplitude and phase demodulation purpose, is very often referred to as complex 

demodulation. Thus'the modulated signal is multiplied by 2e~ Jm * 1 , giving 

(43) 

Then the procedure to obtain demodulated signals fall onto the following steps 39 : 
( a ) determine carrier frequency from the spectrum and bandwidth Aa> occupied by side- 
bands, (b) multiply time series by 2e~ Jeoot , (c) apply a digital low-pass filter and (d) compute 
the amplitude and phase of the complex signal obtained from 3. 

More details about this method can be found in 40 . Either in communication systems or 
in signal processing complex demodulation has received relatively little attention. The main 
difficulty is related to the carrier frequency determination. The carrier frequency asynchronous 
destroys harmonic relations in demodulated spectra. 

4.3 Hilbert Transform Approach 

The analytic signal used in section 3.2.4 to define the envelope function is the basis 
of this approach. Going back to Eq, (22) one can write 

(44) 
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where a(t) is the envelope described in section 2.2.4 and $(t) is the total instantaneous 
phase of the signal being, 

e(t)= arctan =£0^ + <f>{t) 
x(t) 

(45) 

Taking into account the method of the digital amplitude demodulation described in 
section 3.2.4, the phase demodulation can be realized. Both operations can be realized 

37 21 3d 9 

simultaneously (Fig.4). This algorithm is well established and often used in practice ■ 
The method requires removing the carrier frequency part from the total instantaneous phase. 
This can be done by translating the spectrum components to the negative frequencies, so that 

« the component originally at ^ will be moved to zero. Since in practice the carrier frequency 

is determined with some error instead of ^(V) one receives ^(0+4~T " U$in9 th * 

Fourier transform and multiplication theorem of the Fourier transform 28 , one will receive in the 
frequency domain the spectrum of the instantaneous phase 
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(46) 



which gives an unwanted component in the low frequencies. There exist three methods to 
improve the results. First, the differentiation operation of the total instantaneous phase d(t) 
can be performed which gives us the instantaneous frequency /(*) - This operation is not 
effective at all. The other two ways of removing low^frequency trends are 33 least-squares 
polynomial trend removal and high-pass filtering. Since the phase modulation processes are 
low-frequency type, using high-pass filtering, one should be very careful in order not to 
remove the required data. In practice the polynomial trend removal is preferable 33 . 



4.4. Zero-crossing Detection 
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The zero-crossing frequency can be described on the basis of Rice's frequency 
defined as 19 



\ 



(47) 

From the physical point of view this value indicates the central frequency on which 
the whole power is concentrated. For a narrow-banded Gaussian process Rice's frequency is 
equal to the expected rate of zero-crossings with positive slope, which can be written 

E[nJq)]= Km — 
(48) 

where N is a number of "positive" zero-crossing values and At is the interval of time. In 
practice Eq. (48) can be replaced by 

(49) 

Thus practically, it resolves into calculations of the interval of time for which the 
change of the phase has a 2n value, which gives us an estimated instantaneous frequency. 
This method used in VA diagnostics 21 requires a great number of samples in one carrier 
period and does not give sufficient accuracy. However, it does not require the carrier 
frequency determination. 



4.5. Wigner-Ville Distribution 
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The Wigner distribution (WD) can be derived by generalising the relationship between 
the power spectrum and the autocorrelation function for non-stationary, time-variant 
processes. The physical interpretation of the generalised power spectrum f{z y f) is that it 
represents the instantaneous power density spectrum. This leads to the WD defined as 41 , 



(50) 

It can be shown that the first derivative of an arbitrary signal's phase reflects the 

=3 mean instantaneous frequency. The WVD can also be used to extract the instantaneous 

O 

O frequency: 

5*5 s 

CI Let x(f)**A(f)e 1 *® where A(t) and ^(/), the magnitude and phase respectively are real 

valued functions, the first derivative of the phase is given by 41 , 



27V J 



(51) 



This simply means that at time t the mean instantaneous frequency of the signal is 
equal to the* mean instantaneous frequency of the WVD. The disadvantage is that at any time 
instant t there is more than one frequency component present, i.e. the instantaneous 
frequency is not a single value, signal energy spreads with respect to the mean instantaneous 
frequencies. 



4.6. Wavelet Transform: Zero-crossings Detection 
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The wavelet transform is used to decompose a signal x(t) into wavelet coefficients 
l{Vy,x)(a 9 b) using the basis of wavelet functions y/ a ^ if). This can be expressed as, 

+00 

(^x)M)- \x{t)Y Qib {t)dt 

(52) 

where is the complex conjugate of y/ (.). There are many functions that can be used 

as the wavelet basis functions, an example is the Morlet wavelet 42 . 



The square of the modulus of the wavelet transform can be interpreted as an energy 
density over the (a,b) time-scale plane. The energy of a signal is mainly concentrated on the 
time-scale plane around the so called ridge of the wavelet transform. 



L It can be shown that 1 *' 15 , 

\* (53) 

where 0 X and ^ denote the instantaneous phases of the signal and the wavelet transform 

respectively. The ridge of the wavelet transfomi is directly related to the instantaneous 
frequency of the signal. The ridge is defined as, 

<V(o) 



r(bh m 



(54) 



This is used to obtain the instantaneous frequency directly from the ridge of the 
wavelet transform. There exist two algorithms based on the amplitude and phase of the 
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transform 43 . Application examples include damping estimation procedures 44 ^* 5 and 
identification of nonlinear systems 44 . 

5. EXAMPLES AND APPLICATIONS 



Two Simple examples are shown here to illustrate the application of the instantaneous 
characteristics. The first example shows the power of the Hilbert transform for damage 
detection in gearboxes. The second example from the area of system identification shows 
that often wavelet analysis is better when the instantaneous frequency rather than phase is 

required in calculations. 



5.1 Fault Detection in Gearboxes 

Fault characteristics in gearboxes appear in the modulated non-stationary form of 
impacts, which are exhibited in the envelope and instantaneous phase of the vibration data. 
The first example involves the analysis of vibration data from a spur gear. The data came 
from a simple test rig comprising an input gear with 24 teeth driven by an electric motor and 
meshing with 16 teeth of a pinion 15 The rotational frequencies of the wheel and the pinion are 
25 and 37.5 37.5 Hz , respectively. This results in the meshing frequency of 600 Hz. The 
analysed fault is the loss of part of the tooth due to breakage at a point of the working tip. 
Simply, 1 mm of the facewidth was completely removed. Fig. 2 shows an example of power 

spectra representing normal and damage conditions. The damage condition spectrum 
displays a clear pattern of sidebands around the meshing harmonics. Examples of signal 
demodulation results are given in Fig, 3 and Fig. 4. Time domain averaged signal for damage 
and normal meshing vibration are given in Fig, 3a and Fig. 4a, respectively. Local tooth fault 
exhibits an impulse in the envelope function, as shown in Fig. 3b. Also, the change of phase 
can be observed in the instantaneous phase function in Fig. 3c. In contrast the instantaneous 
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characteristics for the normal meshing vibration remain smooth and do not display any 
disturbances, as shown in Fig. 4b and Fig. 4c. 



5,2. Identification of Nonlinear Systems 
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Fig. 5 shows the impulse response function and its spectrum for the vibration 
response from a single degree of freedom rig nonlinear test rig with a cubic stiffness 
characteristic 46 . 

The amplitude of the wavelet transform is given in Fig. 6. The ridge of the wavelet 
transform was computed to obtain the envelope and instantaneous frequency functions for 
the vibration response. Finally, the backbone curve of the analysed system was constructed 
to give the result shown by the solid line in Fig. 7. This characteristic clearly displays cubic 
stiffness nonlinearity. The dashed line in Fig. 7 gives the backbone curve calculated using the 
Hilbert transform approach. Although the oscillations, which are due to differentiation 
procedure, can be removed, the results clearly show the advantage of the wavelet analysis 
over the classical Hilbert transform approach. 

6. CONCLUSIONS 



The envelope and instantaneous frequency/phase functions, are often used for 
damage detection in rotating machinery, for identification of nonlinear systems and in 
deterministic approach to stationary. 



A number of algorithms for the amplitude and frequency/phase demodulation have 
been briefly described. Finally we are in position to make a comparative analysis. The points 

to be compared are: 

• initial filtering - affects the narrow-band characteristics of the signal; 
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• determination of the carrier frequency - important since the carrier frequency is not 
always known or determined with sufficient accuracy; 

• integration or differentiation - both operations require a sufficiently high sampling 
frequency in order to obtain the required accuracy; 

• initial amplitude limitation - difficult in numerical implementation. 

The comparative results are presented in Table 1. The Hilbert transform and energy 
distribution algorithms seem to be very attractive. The Hilbert transform method is well 
established and easy to implement using the classical FFT algorithm. However, it is only valid 

a .. 

for single component signals. Also, it requires numerical differentiation for the instantaneous 
fy frequency, which is not an easy task. New developments in the area of wavelet analysis can 

*y overtake these drawbacks but often lead to expensive computations. Altogether, it appears 

^ that the application and implementation very much depend on the problem under 

s consideration. 

j*j The envelope and instantaneous frequency/phase functions, are often used for 

O damage detection in rotating machinery, for identification of nonlinear systems and in 

□ 

|== deterministic approach to stationarity. The paper also shows that the envelope and 

instantaneous frequency functions are an important link between classical Fourier approach 
based methods and time-variant analysis. 
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Figure 1 . 

Amplitude and phase modulation signal and its spectrum characteristics. 
Figure 2. 

Power spectra for spur gear vibration data: (a) normal condition (b) damaged tooth. 
Figure 3, 

Signal demodulation characteristics for spur gear vibration data representing normal 
^ condition: (a) time domain averaged signal (b) envelope (c) instantaneous phase. 
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Figure 4. 

Signal demodulation characteristics for spur gear vibration data representing damaged tooth: 
(a) time domain averaged signal (b) envelope (c) instantaneous phase. 



J Figure 5. 

o 

Impulse response function (a) and its spectrum for the analysed nonlinear system. 
Figure 6. 

Wavelet transform of the impulse response function. 
Figure 7. 

Backbone curve of the nonlinear system estimated using wavelet analysis ( ) and the 

Hiibert transform ( — ). 
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Table 1. 

A summary of signal demodulation algorithms. 
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Carrier 
frequency 


Integration 
required? 


Differentiation 
required ? 


Filtering 
required ? 


Comments 






required ? 










Complex 


Yes 


No 


No 


Yes 






Demodulation 














Alqorithm 














Local maxima 
distribution 


No 


Yes 


No 


No 


Problem with time 
of averaging 




Energy 

Distribution 
envelope 


No 


No 


Yes 


No 


Differentiation can 
be replaced by 
integration 
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Hilbert 
transform 
algorithm 


No 


No 


No 


No 


Well established 

uses the FFT 
aigoniniTi oui veiiiu 

only for single 
component analysis 
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Yes 


Yes 


No 


No 
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Wigner-Ville 
Distribution 


No 


No 


No 


Yes 


very sensitive to 
noise in the data 
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Wavelet 
Transform 


No 


No 


No 


No 


suitable for 
multicomponent 

signals but 
computationally 

expensive 
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